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ABSTRACT 

We consider accretion disks consisting of counter-rotating gaseous components with 
an intervening shear layer. Configurations of this type may arise from the accretion of 
newly supplied counter-rotating gas onto an existing co-rotating gas disk. For simplicity 
we consider the case where the gas well above the disk midplane is rotating with angular 
rate +Q and that well below has the same properties but is rotating with rate —Q. Using 
the Shakura-Sunyaev alpha turbulence model, we find self-similar solutions where a thin 
(relative to the full disk thickness) equatorial layer accretes very rapidly, essentially at 
free-fall speed. As a result the accretion speed is much larger than it would be for 
an alpha disk rotating in one direction. Counter-rotating accretion disks may be a 
transient stage in the formation of counter-rotating galaxies and in the accretion of 
matter onto compact objects. 



Subject headings: accretion disks - counter-rotating galaxies 
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1. INTRODUCTION 

The widely considered models of accretion 
disks have gas rotating in one direction with a 
turbulent viscosity acting to transport angular 
momentum outward (Shakura & Sunyaev 1973). 
However, recent observations point to more com- 
plicated disk structures in both active galactic 
nuclei and on a galactic scale. Warped and tilted 
nuclear accretion disks have been detected for ex- 
ample in NGC 4753 ( ^teinman-Cameron et ah 



1992| ). Recent high spectral resolution studies of 



normal galaxies has revealed counter rotating gas 
( piri, Bettoni, &: Galleta 1995| ) and/or stars in 
many galaxies of all morphological types - ellipti- 
cals, spirals, and irregulars (see reviews by Rubin 
1994 and Galletta 199 6). NGC 4826 praun et 



al. 19941 ) and IC 1459 ( [Franx fc Illingworth 19881) 
are examples of galaxies with central counter- 
rotating gas and stellar disks, respectively. In el- 
liptical galaxies, the counter-rotating component 
is usually in the nuclear core and may result from 
merging of galaxies with opposite angular mo- 
mentum. Newly supplied gas with misaligned an- 
gular momentum in the nuclear region of a galaxy 
may have important consequences for nuclear ac- 
tivity if there is a rotating massive black hole at 
the galaxy's center (Scheuer 1992). In contrast, 
in a number of spirals and SO galaxies, counter- 
rotating disks of stars and/or gas have been 
found to co-exist with the primary disk out to 
large distances (10 — 20 kpc), with the first exam- 
ple, NGC 4550, discovered by Rubin, Graham, 
and Kenney (1992). It is not likely that the large 
scale counter-rotating disks result from mergers 
of flat galaxies with opposite angular momenta 
because of the large vertical thickening observed 



in simulation studies of such mergers ( Barnes 
1994). Thakar and Ryden (1996) discuss dif- 



ferent possibilities, (a) that the counter-rotating 
matter comes from the merger of an oppositely 
rotating gas rich dwarf galaxy with an existing 
spiral, and (b) that the accretion of gas occurs 
over the lifetime of the galaxy with the more re- 
cently accreted gas counter-rotating. Subsequent 



star formation in the counter-rotating gas may 
then lead to counter-rotating stars. The two- 
stream instability between counter-rotating gas 
and co-rotating stars may enhance the rate of 
gas accretion ( [Lovelace, Jore, fc Haynes 1996| ). 

An important open problem is how counter- 
rotating gas disks form and what their structures 
are on galactic scales and on the scale of disks in 
active galactic nuclei. Here, we investigate accre- 
tion disks consisting of counter-rotating gaseous 
components with gas at large z rotating with 
angular rate +r2(r) and gas at large negative z 
rotating at rate — 0(r). The interface between 
the components at 2;~0 constitutes a supersonic 
shear layer and is shown in Figure |l|. 




Fig. 1. — Structure of two apposed, counter- 
rotating accretion disks and the midplane bound- 
ary layer. The inset shows the three dimensional 
view of the velocity field for the n = 1/2 case 
shown in Figure 2. The velocity variation is anal- 
ogous to that in the Ekman layer of a rotating 
fluid (such as the ocean) where the Coriolis force 
balances the viscous force ( Batchelor 1967| ). 



Configurations of this type may exist in astro- 
physical settings such as halo-disk interactions 
and the accretion of newly supplied counter- 
rotating gas onto an existing co-rotating disk. It 
might at first be supposed that powerful Kelvin- 
Helmholtz instabilities heat the gas to escape 
speed and rapidly destroy the assumed config- 



2 



uration. However, supersonic shear layers exist 
and exhibit gross stabihty in stehar and extra- 



galactic jets ( Hardee, Cooper, &: Stone 1994 ). 
In the counter-rotating disk, matter approach- 
ing the equatorial plane from above and below 
has angular momenta of opposite signs with the 
result that there is angular momentum annihi- 
lation at z = 0, the matter loses its centrifugal 
support and accretes at essentially free-fall speed. 
On the other hand, accretion disks rotating in 
one direction are modeled assuming a turbulent 
viscosity which is crucial for the outward trans- 



port of angular momentum (Shakura & Sunyaev 
19'!^). The counter-rotating disks can also be 
expected to be turbulent owing in part to the 
Kelvin-Helmholtz instability, and turbulent vis- 
cosity can transport angular momentum outward 
in the large \z\ regions of the disk. 

2. THEORY AND SELF-SIMILAR SO- 
LUTIONS 

For stationary, axi-symmetric disk-like flows 
in cylindrical coordinates (r, cf), z) with 

=[vr{r,z),v^{r,z),Vz{r,z)], (1) 
the momentum and continuity equations are 
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Here, p(r, z) is the gas density, p{r, z) is the pres- 
sure, g = — V$ is the gravitational acceleration 



with $(r, z) the potential. For simplicity, terms 
involving the gradient of the dynamic viscosity 
{pu) have been neglected. Because the micro- 
scopic viscosity is negligible for the conditions 
considered, we assume v arises from small scale 
turbulence and can be approximated by the "al- 
pha" prescription of Shakura and Sunyaev (1973) 
as = aCgH, where H is the full half-thickness 
of the disk which is assumed thin (H <^ r), Cs 
is the isothermal sound speed, and a is a di- 
mensionless constant much less than unity. Since 
the turbulent motions may not be isotropic, we 
allow two different tensor components of u (z^y 
for shear in the f — (p direction and for shear 
in the vertical direction) with possible different 
r— dependences. We also assume the second vis- 
cosity coefficient u' 

J I - J i •'J LJ 

Here, v = 

if J = z, and v = z/y otherwise, yielding equations 

We seek self-similar solutions to equations (||)- 
|5|) of the form 

Vr{r,z) = -Ur{C,)Vc{r), 



|i/ so that the stress tensor 
^ tj — pviVj +p5ij + T^j with the viscous contribu- 
tion T^^- = —^{dvi/dxj + dvj/dxi 
z, and v 



Vz[r,z 



v^{r,z) = u^{C)Vc{r), 
h{r) 



Uz{()V,{r) 



(6) 



p{r,z) = po{r)Z 



H(r] 



where C = z/h{r) is the dimensionless vertical 
distance in the disk with h{r) a length scale iden- 

1 /2 

tified subsequently, and Vc(r) = [r{d^/dr)\z=o] 
is the circular velocity of the gas. The potential 
is due in general to disk and halo matter and 
a central object. Here and subsequently we ne- 
glect the pressure force in the radial equation 
of motion. We consider cases where Vc{r) = 
{GM/roy/'^{ro/r)"', with rg a constant length 
scale. The value n = 1/2 corresponds to a Ke- 
pler ian disk around an object of mass M, and 
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n = to a flat rotation curve applicable to flat 
galaxies. 

Substitution of equation @ into equations @ 
and (0) gives 



/I I 2 2/-2\ // 2 I 2 1 I A / 

[1 + e 7] Q )Uj. = + nu^ — 1 + r]QurU^ 

-u^u'^ - e^Cvir] + 2nX + e^{l - r?)ur, (7) 

(1 + e^r]'^Q^)u'^ = (n - l)urU^ - u^u'^ 
+vCuru'^ - e^Cviv + 2n)u; + e^{l - r?)u^, (8) 

where a prime denotes a derivative with re- 
spect to C- We have chosen the scaling such 
that v^r/{Vch?) = 1 and = i>\\/rVc-, and de- 
flned rj = (r /h)(dh/dr). We assume that the 
disk is not strongly self-gravitating and that it 
is vertically isothermal so that its overall half- 
thickness is H {cs/Vc)r which is the same 
as for a disk rotating in one direction. From 
the assumed scaling, /i^ = i'±r/Vc, and an as- 
sumed alpha viscosity prescription z^y = aCgH 
and u± = acsh (|Regev 19831) , 6 = h/H a '^l. 
As a result, @ simplifles to (dzp) = —pdz^. 
which determines the isothermal density proflle 
Z{() = exp(— Prom the assumed scaling 
Po{r) oc r~'^~^/^, equation 



P + n 



2 a2 



26' C 



r dH' 



(9) 



determines Vz ■ For consistency of the self-similar 
solutions, e must be independent of r, implying 
rj = 1 and [r / H){dH / dr) = 1 if the viscosities 
scale similarly, v±_ ^ 

Equations (0)-(|) constitute a closed system 
for the dimensionless functions Ur{C,),u^{C,), and 
Uz{C)- Applying equations to the flow 

suggested in Figure 1 |l|, we infer that Ur(C) is an 
even function of Q, while u^{C,) and UziQ) are odd 
functions. For |^| ^ 1, we impose Ur ~ (1 -|-n)e^ 



and ~ ±\/l — {n + l)^e^ so as to have u'l — > 
and u'^ 0. These limits correspond to an a 
disk rotating in one direction far away from the 
C = midplane. For n = 1/2, ii,.(C » 1) ~ 3e/2 
and Vr{C ^ 1) = 3z/||/2r, agreeing with the ra- 
dial velocity of a single thin Keplerian disk with 
viscosity = e^rVc- 

A steady counter-rotating disk may result from 
gas continuously supplied at large r. Similarly, 
steady state solutions where gas is entrained and 
delivered vertically onto the faces of the disk may 
exist. We assume that Vzp{r, z) ^ as \z\ — > oo. 
Considering only this subset of solutions, we de- 
mand the accretion rate 



M = 2tt dzrp{r, z)\vr{r, z)\ ^ 2TTrp(){r) 

J —oo 

xV,{r)h{r) f"^\cur{C)e-'"'^\ (10) 

J-H/h 

be a constant. A sufficient condition for this is 
fi = 7] — n + 1/2 and 6 = h/H is r— independent. 
For an alpha disk rotating in one direction, Mss = 
2TTrT,\vr\sSi with the accretion speed jfrlss ~ 
aCgH/r <C c^, where S = / dzp is the disk's 
surface mass density ( Shakura &: Sunyaev 1973 | ) . 
In contrast, for a counter-rotating disk equation 
( [lOD gives McR = 2-KrTi\vr\cr 1 with accretion 
speed (averaged over z) \vr\cR ~ {h/H)Vc. For 
the same a, we have |vr|ci? ^ brls^: the accre- 
tion speed of the counter-rotating disk is much 
larger than that of the disk rotating in one direc- 
tion. 

A further constraint on solutions of equations 
(|2|)-(^ is obtained by considering the energy dis- 
sipation and radiation of the disk. The viscous 
dissipation per unit area of the disk is 



D{r) = 

<^ 

where Q 



H 



H 



+ 



dz p V 



(dv. 



VA r. 



\ dr _ 

(11) 

For a disk rotating in one 



dvz ^ dvr_y ^ 
dr dz 
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direction, the dominant contribution to D{r) is 
from the dQ,/dr term and this gives Dssif) ~ 
Sz.||(n + l)2(K/r)2 = Mss{Vjry{n + l)V(37r). 
In contrast, for a counter-rotating disk, the dom- 
inant contributions are from the terms dv^/dz 
and dv^jdz with the result 

~ Mcfl(Fe/r)V(37r) ~ {TjiifDssir). (12) 

The dissipated energy is radiated from the faces 
of the disk if it is optically thick in the z-direction 
and thus D = 4acT^/ (3kS), where a is the radi- 
ation constant, c the speed of light, k the opacity, 
and T is the internal disk temperature. For the 
general scaling k oc p°'T^, the counter-rotating 
disk thus has T oc r^^ with ^ = (3 + 2a -|- n)/(4 — 
h + a/2) and h oc qc r"^ . For example, 

for Kramer's opacity a = 1, h = —3.5, we have 
^ = (5 -|- n)/8 and 77 = (11 -|- 15n)/16, whereas 
for electrons scattering opacity a = 0, 6 = 0, and 
^ = (3-|-n)/4 and r/ = (ll-|-9n)/8. In general, the 
strict requirement r] = 1 for our self-similar solu- 
tions will not be satisfied, but this is not expected 
to be important if e <C 1 (see Figure ||)[|. The ap- 
parent surface temperature of the disk Teff{r) is 
given by D{r) = 2aT^jj, where a = ac/A is the 
Stefan-Boltzmann constant. The dependence for 
the counter-rotating disk, Tf.ff oc jg 
same as for an alpha disk rotating in one direc- 
tion. Similarly, the spectrum, F^^, obtained by 
integrating the Planck function B^{Tgff{r)) over 
the surface area of the disk, F^^ oc u}^^"'~^^^^"'~^^\ 
is the same as for a disk rotating in one direc- 
tion. For an optically thin disk, D = 2HA, where 
A{p, T) is the emissivity of the disk matter. 

In the limit e = 5 = and (3 = rj — n + 1/2, 
equation (^) gives Uz = TjQur-, which reduces 

^ Because the viscosity scaling ~ z/y implies iq — 1, 
values of 7^ 1 are strictly consistent with our scalings and 
constant M only if the turbulent viscosity is anisotropic 
such that shear perpendicular to the plane is governed by 
— e^rVc ~ but shear in the plane is determined 

by u± = h?Vc/r ^ 



equations (j^ and (^) to u'l = v?^ + nu^ — 1 and 
u'^ = {n — l)urUtj,, respectively. These have the 
integral 

K)' - ^^«)' + (1 - ^)^rn^ 

n(l — n) n , , , ^ 

H — ^— -ui - (1 - n)ur = const. (13) 

Evaluating (|l^) at ^ = and C, = oo gives 
[u'^(0)]2 = (1 - n)ur{Q) - n(l - n)u3(0)/3. We 
utilize this result by first finding the single inde- 
pendent initial condition Ur{Q) for (5 = e = 0, 
then allow 5 , e 7^ 0, and perturbatively find 
the proper values of ^,.(0) and ti^(O). Three 
of five initial conditions are fixed by symmetry: 
ti^(O) = u^{Q) = Uz(S^) = 0. We search for so- 
lutions by tuning tir(O) and m^(0), regarding e 
and 6 as fixed. In principle, the required val- 
ues u,.(oo) and u^{<X)) are satisfied for a range of 
5. The additional parameter (5 ^ e can be cho- 
sen according to the behavior of Uz{oo) desired. 
Thus, a class of solutions exist located on curved 
segments in [^^(0), u^(0), 5] space. 

Figure ^ shows solutions to equations (|^)-(|9|) 
for e = 0.01 and 5 = 0.1, corresponding to 
- 1/10 and /i/r ~ 1/300. Both Keplerian 
(ra = l/2,/3 = 1) and galactic (n = 0,/? = 3/2) 
solutions are shown. The solutions are sensitive 
to the precise values of tir(O) , u^(0) , 5, and e. 
Note that for 3 ^ |C| ^ 7, the solutions have a 
region with ti^ < indicating spiraling outflows 
on both sides of the midplane. The vertical ve- 
locities Vz/Vc ^ Ur or are indistinguishable 
from the x-axis on the scale of Fig. ^. For large 
C, the density fall off gives pvz 0. 

3. DISCUSSION 

Counter-rotating gas supplied to the outer 
part of an existing co-rotating gas disk will in- 
crease the mass accretion rate. Comparing ac- 
cretion rates of a standard a— disk (Mss) with 
that of a counter-rotating Keplerian disk (Mcr) 
of the same S and M, we find from equation 
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1.4 H 




Helmholtz instabilities are likely to be the most 



important ( Ray fc Ershkovich 1985 and Choud 



Fig. 2. — Solutions to Equations (||), and 
(^) for e = 0.01 and 6 = 0.1 as a function of 
C, = z/h{r). Solid lines correspond to a Ke- 
plerian disk (n = 1/2, fi = 1) with ^^(0) = 
1.10007597 . . . and n'^(O) = 0.66260669 .... The 
dashed curves correspond to n = (flat rota- 
tion curve potential) and /3 = 3/2 and have 
initial values ^^(0) = 1.26728681248... and 
u'^{0) = 1.1256582571 .... These values are close 
to those of the 6 = e = "zero-order" solutions: 
Ur{0) = K^(0)]2 = 1.100112673... (n = 1/2) 
and Ur{0) = [u'^iO)]^ = 1.2672806222 ... (n = 0). 
All solutions reach their asymptotic values for 
^ oo and the vertical velocity Vz/Vc ~ on 
this scale. Note that the midplane infall velocity 
for the Keplerian case is shear-reduced by a fac- 
tor ~ 0.77 from the free-fah speed (2GM/r)i/2_ 

(|10D Mcr/Mss ~ 1 + 0.58(5/e2). For the values 
of Figure 2 this ratio is ~ 580. 

Thermal and/or dynamical instabilities may 
destroy the counter-rotating accretion flows de- 
scribed above. The relative importance of ther- 
mal and dynamical instabilities can be estimated 
by comparing the thermal dissipative time scale, 
TQ « J:cyDcR{r) ~ {r/h){H/Vc){cs/Vcf, with 
dynamical and viscous time scales, Tz = h/cs, t± = 
r/Vc, and Ti, = jv^. We find tq ~ e^b~'^Tz ~ 
e^^-^/Vx ~ e'bT^. For e = 0.01 and b = 0.1, the 
thermal dissipation time tq <^ Tz < t± <^ t^- In 
contrast to the standard a— disk, thermal equilib- 
rium is maintained on time scales of possible flow 
instabilities of the inner shear layer. The Kelvin- 



hury fc Lovelace 1984| ). If the counter-rotating 
disk is treated vortex sheet, then a local 

stability analysis indicates unstable warping for 
wave numbers \k(f,/kr\ < V2{cs/^r) <C 1. 

Accretion of counter rotating gas by an ex- 
isting co-rotating gas disk may be a transient 
stage in the formation of counter-rotating galax- 
ies and in the accretion of matter onto rotating 
black holes in active galactic nuclei. We find that 
newly supplied counter-rotating gas drags inward 
the old co-rotating gas with an equal mass of old 
and new gas accreting rapidly. Thus the old co- 
rotating gas may be entirely "used up" (dragged 
to the center of a galaxy or into a black hole) if 
the mass of newly supplied gas exceeds that of 
the old gas disk. Accretion onto the faces of an 
existing thin disk may not have the symmetry 
shown in Figures Q and ^. 



® 



Fig. 3. — Schematic drawing of accretion of 
newly supplied counter-rotating gas (©) induced 
by viscous interaction with an existing disk of 
co-rotating gas (©). 

There may instead be two layers of rapid radial 
inflow bounding the existing gas disk near the 
midplane as sketched in Figure ^. Solutions for 
this configuration can be composed from those of 
Figure | if /i<C 
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